Electron transport and current fluctuations in short coherent conductors 
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Employing a real time effective action formalism we analyze electron transport and current fluc- 
tuations in comparatively short coherent conductors in the presence of electron-electron interac- 
tions. We demonstrate that, while Coulomb interaction tends to suppress electron transport, it may 
strongly enhance shot noise in scatterers with highly transparent conducting channels. This effect of 
excess noise is governed by the Coulomb gap observed in the current- volt age characteristics of such 
scatterers. We also analyze the frequency dispersion of higher current cumulants and emphasize a 
direct relation between electron-electron interaction effects and current fluctuations in disordered 
mesoscopic conductors. 
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I. INTRODUCTION 

Recently considerable progress has been reached in un- 
derstanding of an interplay between scattering effects and 
electron-electron interactions in low-dimensional disor- 
dered conductors. In particular, a profound relationi 2 *^ 
between full counting statistics 5 (FCS) and electron- 
electron interaction effects in mesoscopic conductors has 
been discovered and explored. This relation follows di- 
rectly from the effective action of the system which 
can be conveniently derived by combining the scatter- 
ing matrix technique with the path integral description 
of interactions 1 - 2 - 6 . 

An important simplification of the formalism amounts 
to neglecting the energy independence of scattering ma- 
trix. This approximation is applicable to rather short 
coherent scatterers with Thouless energy exceeding all 
other relevant energy scales. Besides that, all inelastic 
processes are assumed to occur in large reservoirs but not 
inside the scatterer. For this model an exact expression 
for the effective action was derived 3 *^ and studied under 
various approximations, such as regular expansionist of 
the action in powers of the "quantum" part of the fluc- 
tuating field describing interactions as well as renormal- 
ization group (RG) analysis^. Both approaches are jus- 
tified in the metallic regime, i.e. provided the effective 
system conductance is much larger than the quantum 
unit e 2 jh. Further generalizations of the model allowed 
to consider spatially extended metallic conductors and to 
describe inelastic processes inside the systeirA 

The main goal of the present paper is to investigate 
the problem employing a different set of approximations 
in order to go beyond the regimes already studied in the 
literature. Having in mind the relation between FCS and 
interaction effects in mesoscopic conductors, it would be 
useful to develop a straightforward and handy formal- 
ism enabling one to analyze the frequency dependence 
of higher order current correlation functions. In other 
words, the task at hand is to generalize the FCS-type of 
approach 5 (which is valid only in the low frequency limit) 
to arbitrary frequencies. In addition, it is desirable to ex- 
tent the description of interaction effects in mesoscopic 



coherent conductors beyond the metallic regime. This 
step would allow, for instance, for a more detailed analy- 
sis of electron-electron interactions effects in conductors 
with few conducting channels, such as, e.g., single- wall 
carbon nanotubes and organic molecules. 

Here we will demonstrate that the above goals can be 
accomplished assuming transmissions T n of conducting 
channels to be either sufficiently small or, on the con- 
trary, sufficiently close to unity. In these situations sig- 
nificant simplifications of the exact effective action can be 
worked out making the whole formalism more convenient 
for practical calculations. 

Perhaps the most striking result of our analysis con- 
cerns the effect of electron-electron interactions on the 
shot noise in conductors with several highly transmitting 
channels. We will demonstrate that in this case Coulomb 
interaction yields strong enhancement of the shot noise. 
Specifically, at sufficiently small voltages V the low fre- 
quency noise power spectrum 5/(0) in weakly reflecting 
scatterers takes the form 



5/(0) <x {Vl 1 - 



(1) 



where g is the dimcnsionless conductance of a scat- 
terer. Provided g is not much bigger than 2 the result 
strongly exceeds the non-interacting dependence^ 
5/(0) oc \V\. In the limit of large voltages the noise 
spectrum is offset by the value of the Coulomb gap, 



5/(0) oc V + e/2C, 



(2) 



where C is the effective capacitance of the scatterer. We 
also note that interaction-induced excess noise in weakly 
reflecting scatterers is observed only at sufficiently low 
frequencies w < eV, while at larger uj Coulomb interac- 
tion - on the contrary - suppresses the shot noise. 

The structure of the paper is as follows. In Sec. II we 
will briefly summarize the main steps of our derivation of 
the exact effective action for the conductor described by 
an arbitrary energy independent scattering matrix. This 
general expression will be analyzed in Sec. Ill in various 
limits. The limiting forms of the effective action will then 
be used to describe electron transport and current fluc- 
tuations in coherent conductors. In Sec. IV we will use 
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them to obtain the frequency and voltage dependence of 
current cumulants. In Sec. V we will consider a scatterer 
shunted by a linear Ohmic conductor and, assuming that 
the system conductance is sufficiently large, discuss per- 
turbative logarithmic interaction corrections to the I—V 
curve as well as their RG analysis. In Sec. VI and VII 
we will go beyond the metallic limit and study the effect 
of electron-electron interactions respectively on the I — V 
curve and on the shot noise of highly transmitting con- 
ductors. Further extensions of our results for the case of 
quantum dots will be considered in Sec. VIII. Some for- 
mal manipulations and further details of our derivations 
are presented in Appendices A, B and C. 



II. EFFECTIVE ACTION 

For our derivation we will make use of the Keldysh path 
integral formalism describing the system of interacting 
electrons by the action which depends on the Grassmann 
electron fields *£.' 2 (i, r), #^ 2t (t, r). The labels 1 and 2 
correspond to the forward and the backward branches 
of the Keldysh contour and a is the spin index. Per- 
forming the standard Hubbard-Stratonovich decoupling 
of the Coulomb interaction terms in the action, one intro- 
duces the two fluctuating electric potentials V\ (t, r) and 
V 2 (t, r). Afterwards the electron fields are integrated out 
and one arrives at a formally exact effective action 

iS = 2THnG- + i f dt> f d*r ™* = ™* . (3) 
Jo J 8tt 

Here we assumed the spin degeneracy and defined 











iTt+h + ^-U{r) + eV 2 



,(4) 



where U(r) is the external static potential. 

The action can also be expressed in a different way, 
which is convenient for generating correlation functions 
of the current: 



iS = In Tr 



f dt'H x (t') « ^- i f* dt'H 2 (t') 



8?r 



(5) 



with the trace taken over the fermionic variables. Here 
Po is the initial N— particle density matrix of electrons, 

H h2 (t) = -^- + U(r)-eV 1 , 2 (t,r) (6) 

are the effective Hamiltonians on the forward and back- 
ward parts of the Keldysh contour, T and T are respec- 
tively the forward and backward time ordering operators. 



One can also show (see Appendix A for details) that 
the action @ can be expressed in the following equivalent 
formS 



iS = 2trm{l + (W 2 (0,t)W 1 (t,0)-l)A)} 



i J dt' i <!'•>■ 



(Wi) 2 - (W 2 
8tt 



where 



U 



+ n ^ ^f f 2 dt'(-Sl~ l ,+U(r)+eV 1 , 2 (t\r)) 



(7) 



(8) 



are the evolution operators. Here we distinguish between 
the "full" trace ("Tr", Eq. implying the integration 
over both the coordinates and the time variable, and the 
"reduced" trace ("tr", Eq. (JJJ) which denotes the coor- 
dinate integration only. 

In this paper we will analyze the properties of a short 
coherent conductor placed in-between two bulk reser- 
voirs. Electron transport across such a conductor can 
be described by the scattering matrix: 



S = 



f t' 

i ? 



(9) 



Below we will assume that this matrix does not depend 
on the energy of incoming electrons. This assumption is 
justified provided the dwell time of an electron inside the 
scatterer is shorter than any other relevant time scale. 
In addition, we assume that an external circuit, which 
also includes the leads, is linear and characterized by an 
effective impedance Zs{uj). The fluctuating electric po- 
tential is assumed to slowly depend on the coordinates 
in the reservoirs, but to drop sharply across the scat- 
terer. Denoting this voltage drop as Vi^{t r ), one can 
also introduce the corresponding Keldysh phase fields 
<Pi,2(t) = L dt' eVifi[t') as well as their combinations 
ip + = (ipi + tp 2 )/2 and ip~ = ipi — ip 2 . The total ac- 
tion of our system reads 

iS[tp ± ]=iS s [<p ± } + i db'C^- + iS sc [<p ± ]. (10) 
Jo e 

Here iSs[tp ] describes the external circuit: 



^ Jo 



dt 1 dt 2 <p (ti)a s (tx -t 2 )tp (t 2 ) 



+ 4 / dt 1 dt 2 ^-(t 1 )Z s 1 (t 1 -t 2 )(eV x -^+{t 2 )), {11) 
e Jo 

where Zg 1 (t) = J ^ z s {u) ^ s ^ e res P onse function of the 
external shunt, a s (t) = J % e'^Re ( "^j* ) , and 
V x is the total voltage drop on the system "scatterer + 
external shunt" . 

The second term in Eq. (|10fl originates from the last 
term of Eq. J7J) and describes the energy of the fluctuat- 
ing fields. This term contains the scatterer capacitance 
C. 



3 



Finally, the most interesting for us contribution iS sc 
is given by the first term of Eq. (JJJ with the evolution 
operators U1.2 describing the scattering of electrons in 
the conductor. Following the standard procedure we in- 
troduce the scattering channels (labeled by indices) and 
denote the coordinate and the group velocity in the k— th 
channel respectively as y k and v k . The matrix corre- 
sponding to the operator U\(t, 0) then takes the form 

s(m.- m. - t ) 

U? k (t,0;y 2 , yi ) = e **»yi(t)-m fc yi(o) Vj; / x 

y/v n v k 

[s nk + e{y 2 )6(~ yi )e- m ^-y^ a ^ { ~^\s nk - S nk )] , 



(12) 

where S nk is the matrix element of the S— matrix 
a k = 1 provided the k— th channel is in the left lead and 
a k = otherwise. The operator lA2{Q,t) is constructed 
from Ui(t, 0) by means of the Hermitian conjugation to- 
gether with the replacement (pi — > if2- We point out that 
Eq. (|12|l is similar to the well known relation between the 
Green function and the S— matrix established by Fisher 
and Leeii. 

Having found the expressions for the evolution oper- 
ators, we make use of Eq. Q and arrive at the final 
result: 



iS sc = 2trln <^ 15(x - y) + 6(t - x)9{x) 



>itp (x) 



1) 



2^r't^sin^^ t*£( e -*v (*) 



lifir' sin 
(e 



2 



1) 



Po(y - x)e 







po(y - x)e 



■p" 1 " (yj -y + (g0 



r 



113] 



where 

f dx e iEx 1 r . , iT 
P ° (X) = J 27TT^7t - 2 S{X) 2anh*T X (14) 

is the Fourier transform of the Fermi function. 



III. LIMITING FORMS OF THE ACTION 

The general expression for the effective action (|13[) con- 
tains a great deal of information and can be used in or- 
der to describe a variety of different phenomena. At the 
same time this expression still remains rather complex, 
and further simplifications are highly desirable. These 
simplifications can be achieved in various limiting cases 
to be discussed below in this section. 

Before we turn to concrete calculations let us note that 
several important physical limits of Eq. I|13|) are already 
well known and have been studied in details. For in- 
stance, setting (p~ = const and (p + = eV one reduces the 
action (|13|1 to the FCS cumulant generating function^ for 
a coherent scatterer in the absence of electron-electron in- 
teractions. In this case tp~ plays the role of the so-called 
counting field. 

If one allows for the time-dependent voltage V(t) 
across the scatterer and at the same time keeps ip~ con- 
stant, with the aid of Eq. I|13|) one can describe FCS at 
an ac bias, adiabatic pumping through the conductor and 
related effects 1? . In addition, neglecting the interactions 
but keeping the full time dependence for the phase ip~, 
with the aid of Eq. I|13|l one can fully describe frequency 
dispersion of all current cumulants. For this purpose it 
suffices to perform a regular expansion of the action l|13|) 
in powers of the field ip~ . Frequency dispersion of the 



third current cumulant was analyzed in this way in Ref. 
I 

The same expansion allows to obtain valuable infor- 
mation about the effect of electron-electron interactions 
on transport properties, shot noise as well as higher cu- 
mulants of the current operator. Assuming the dimen- 
sionless conductance of a scatterer g and/or that of the 
external circuit gs are/is large, one can obtain the in- 
teraction correction to the current^ of order l/(g + gs) 
expanding the action (fT~Sf) up to the second order in Lp~ . 
In order to derive the interaction correction to the shot 
noised one should expand the action up to the third order 
in tp~ . 

Finally, keeping the exact non-linear dependence of the 
action on the fields ^ but expanding (|13|l to the first 
order in the transmission matrix fit, one immediately re- 
produces the well known AES effective actioni^i 4 . The 
latter action, being combined with iSs and the ca- 
pacitance term, describes Coulomb blockade effects in 
tunnel junctions embedded in a linear electromagnetic 
environment 1 ^* 1 ^. 

In what follows we will investigate the properties of the 
action l|13|l in some other limiting cases. 



A. Weakly transmitting barriers 

Let us first consider the case of weakly transmitting 
barriers and expand the action in powers of the matrix 
fit. As we have already pointed out, the first order terms 
of this expansion yield the AES effective actioniiiit. Here 
we proceed further and expand the action (|I3|I up to 
(iH) 2 keeping the complete nonlinear dependence on the 
fluctuating phases ip ± . It is easy to see that in order to 
recover all such terms it is necessary to expand the log- 
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arithm in Eq. (|13fl up to the fourth order in the term 
containing the density matrix p . Higher order terms of 
this expansion can be omitted within the required accu- 
racy since they do not contain contributions proportional 
to T n and T„, where T n represent the channel transmis- 



sions defined in a standard manner as the eigenvalues of 
the matrix t't. 

The whole calculation is performed in a straightfor- 
ward manner, although requires some care. The final 
result reads 



iS 



— / dx if (x) sin ip (x) 



(i 



trpi] + 2 -tr[(iH) 2 } sin 2 ^ 



-- / dxdy a(x - y) sin ?-M. sin T-Mi / tr [^t(l - *+*)] cos[^+(x) - <p+(y)] + tr[(Pt) 2 } cos t g ^ 

7T ,/n 2 2 2 



4i 

+— tr[(t+t) 2 ] / cfccfydz 



T 3 sin sin ^ sin 



sinh[7rT(y — a;)] sinh[7rT(a; — z)] sinh[7rT(z — y)] 



x <j sin[<^ + (2/) — iy9 + (a;)] cos - ^ + sin[<y9 + (z) — if + (y)} cos - ^ + sin[(^ + (x) — f + (z)] cos - ^ 



— 16tr[(f^i) 2 ] / dxdydzdw po(y — x)pq(x — w)po(w — z)pq(z — y) 
Jo 

<p~(x) <p~(y) cp-(z) <p~(w) 



x sm ■ 



■ sm ■ 



■ sm ■ 



■ sm ■ 



■ cos 



[if + (x) - ip + (y) + if+(z) - ip + {w)] 



(15) 



where 



a(x)= I ^e" 



> coth — = — - 
2T 



ttT 2 



sinh" nTx 



(16) 



Eq. I|15() represents the complete expression for the effec- 
tive action valid up to the second order in the transmis- 
sions T n . This expression involves no further approxima- 
tions and fully accounts for the non-linear dependence on 
the fluctuating phase fields (p ± . 



B. Refiectionless limit 

Another physically important limit is that of reflection- 
less barriers f = 0. In this limit the action can be 
significantly simplified by means of the exact procedure 
which we outline below. 

To begin with, we notice that in the case of refiec- 
tionless barriers the action lll.'ifl reveals significant sim- 
ilarity to the Luttinger model of ID interacting elec- 
tron gas. The latter model is usually treated by means 
of the bosonization technique. In the case of quantum 
dots and point contacts this method was applied in Refs. 

One can also evaluate an effective action 
for refiectionless barriers directly without employing the 
bosonization technique. Actually an important feature 
of the final result can be guessed even before doing the 
calculation. Indeed, since the RPA approximation is 
knownSS as an exact procedure for the Luttinger model, 
one can expect that in the limit f' — > the action 
should become quadratic in ^ , at least for not very 



large values of the phase. We will demonstrate that this 
is indeed the case if \<p~\ < ir. 

Let us put f' — in Eq. I|13|) . Then the action can be 
split into two parts 

iSjg> = iN ch S [9(t - x)9(x)(e^^ - 1), -<p+] 

+ iN ch So[9(t~x)9(x)(e' 1 ^^ - 1), if+] , (17) 

where A c h is the number of open channels and 



iSo[a, if + ] = 2trln 



1 + a(x)p (y - x)e 



(18) 

The action l|18|) can be evaluated exactly. The details 
of our derivation are summarized in Appendix B. Here 
we only present the final result: 



tS a [a,if+] = 2 / dx[p (0) + 



47T 



ln[l + a(x)] 



+ [ dxdy a(x-y)^ 1 + a ^ 1 + a ^\ (19) 

J 47T 

where po(0) is a large constant, which has a meaning 
of the electron density, and which is later canceled by 
the corresponding contribution of ions. Combining this 
formula with Eq. (|1 T|) we arrive at the main result of this 
subsection 



^sc — 
A ch 



iN, 



cli 



dxW(ip (x))ip + (x) 



I dxdy W (x))a(x-y)W (y)), (20) 
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where W(<p) is the 2tt— periodic function of ip~ , which 
equals to (p~ in the interval — tt < ip~ < tt. Under certain 
conditions the contribution of large phase values, \if~\ > 
7r, can be disregarded, and one can use the Gaussian 
action 



iN, 



ch 



dx(p (x)<p + (x) 



N, 



^ / dxdyip (x)a(x-y)tp (y) 
27r Jo 



(21) 



instead of the exact one <|2(J|I . Within this approximation 
the action for the scatterer with iV c h perfectly transmit- 
ting channels coincides with that for a linear Ohmic re- 
sistor with the conductance e 2 N c \ 1 /ir. It is important to 
emphasize, however, that this approximation ignores the 
action periodicity in the phase space and, hence, becomes 
inadequate as soon as electron charge discreteness turns 
out to be an important effect. 



C. Weakly reflecting barriers 

Let us now assume that the reflection probabilities are 
small, R n = 1 — T n <C 1, but not equal to zero. In this 
case it is convenient to proceed perturbatively expanding 
the action l|13f) in powers of R n or, which is the same, in 
powers of the matrices r', r'K The details of our deriva- 
tion are provided in Appendix C. Expanding the action 
(|13fl to the first order in R n , one finds 



iS Ba = iSj$ + is£>, 



(22) 



where iSs? is defined in Eq. (|20|l and 



dx <p + (x) sin ip (x) 



1Z 

H / dxdy a(x — y) ip (x) shop (y) 



2K 



dxdya(x-y) e ^ + ^)-v + {y)] 



J* dz[T cothTrT(x-z)-T coth-rrT(y-z)]cp- (z) 



X & 

. <p~(x) . <p~(y) 
x sm sin . 



(23) 



As we have already pointed out, this expression is jus- 
tified provided all the channels are weakly reflecting, 
R n < 1. At the same time the parameter 1Z = ^2 n Rn 
needs not to be small should the total number of chan- 
nels in the system be large. We also note that differ- 
ent limiting forms of the action l|15l) and (|20|l - l|23|l can 
be combined if some channels are weakly transmitting 
T n <C 1 while the others have small reflection coefficients 

Rm < 1. 



IV. FREQUENCY DISPERSION OF HIGHER 
CURRENT CUMULANTS 



Let us now make use of the above limiting expressions 
for the effective action in order to describe various prop- 
erties of short coherent conductors. In this section we 
neglect the effects of electron-electron interactions and 
set <p + (t) = eVt, where V is the time-independent bias 
voltage. 

Complete information about all correlation functions 
of the current operator in the zero frequency limit is con- 
tained in the FCS cumulant generating function^ . This 
function, however, becomes insufficient if one is inter- 
ested in the frequency dependence of the current cumu- 
lants. This dependence can in general be recovered only 
from the complete effective action (|13|l . Unfortunately 
the latter appears too complicated to directly proceed 
with the analysis of the n-th current cumulant. The situ- 
ation is significantly simplified in the two complementary 
limits of small and almost perfect channel transmissions, 
respectively T n <§; 1 and 1 — T n <C 1. In these two cases 
one can make use of the limiting forms of the effective ac- 
tion derived in the previous section. The corresponding 
analysis is presented below. 

In the absence of interaction effects one can treat the 
action iS sc as a generating functional for the current cor- 
relation functions 



(/(ii).../(i m )) = 



(ie) m 5 m giSsdeVt.ip ] 



Sip (ti)...5p (t m ) 



(24) 



Here I(tj) are measurable classical currents which com- 
mute with each othesi. On a quantum level, the same 
correlator can be written in terms of the current oper- 
ators I(tj), however the choice of the time ordering be- 
comes important in this case. 

One can also define the m-th current cumulant 
S m (ti, . . . t m ) as an irreducible part of the correlator 12411 : 



S m = (ie) 7 



S m iS sc [eVt,<p- 



5tp (ti)...5ip (t m ) 



(25) 



For classical currents, the first three cumulants are de- 
fined as 

Sx = {I), 

5 2 (ti,t 2 ) = (6I(h)SI(t 2 )) , 

5 3 (ti, h, h) = (SI(h)5I(t 2 )SI(t 3 )) , (26) 

where SI(t) = lit) — (I). The fourth and higher cumu- 
lants take a more complicated form, e.g., 

SiihMMM) = (SI(h)5I(t 2 )SI(t 3 )SI(U)) 

— ^2)^2(^3, £4) — S2\tl,t3}S2\t2, £4) 

-&(ti,t4)5 2 (t 2 ,t3). (27) 

In what follows we will use Eq. (|25|l as a definition 
of the current cumulants also in the quantum case. This 
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definition unambiguously fixes time ordering of the cur- 
rent operators. The corresponding expression for the 
third cumulant (m = 3) in terms of the time-ordered cur- 
rent operators has been specified in Ref. 0. Analogous 
expressions for higher cumulants to > 3 are cumbersome 
and we do not present them here. 

We also define the Fourier transform of the current 
cumulants: 



S,, 



dt\ . . . dt„ 



a iu)it\ + ...+iu) m t n 



2tt 



' S m (tl 



■ in 



(28) 



We begin with the limit of weakly transmitting barriers 
in which case in the lowest order in T n one can evaluate 
the current cumulants making use of the AES effective 
Combining Eqs. (|15fl and l|25[> and dropping 



actio: 



.13.14 



the terms oc T% for the odd (to 
finds 



2^ + 1) cumulants one 



m+l 



S m — 



gv 



2tt 



S(uji + ... +uj m ), 



(29) 



where g — 2 T n is the dimensionless conductance of 
the scatterer. Analogously one can evaluate the even 
(to = 21) current cumulants which read 



£ 

I/,-=zfcl 



eV - 



3 = 1 2 



coth — — 
2T 



\ - VjUj 



3 = 1 



AT 



•(30) 



Here the prime in the sum implies the summation over 
"charge" configurations i/j = ±1 with the odd number of 
positive (negative) "charges" Vj . This result is also valid 
in the lowest order in T n . 

Keeping the terms oc T% in Eq. (fl5|l and repeating 
the same calculation one can evaluate the second order 
corrections to Eqs. I|29[l and (|30(l . The correspond- 
ing expressions turn out to be rather complicated and 
for this reason are omitted here. We note, however, 
that in the limit of low voltages (or high frequencies) 
eV -C T,u>i, . . . , uj m the non-local in time terms in the 
action Ijl5(l do not contribute to the odd current cumu- 
lants. Hence, in this limit the whole analysis gets much 
simpler and the odd cumulants are fully determined by 
the remaining (local in time) part of the action. In other 
words, in the limit eV <C T,u>i,... ,tu m the odd current 
cumulants can be evaluated (up to the terms ~ T^) with 
the aid of the generating function 

ieV 



iSd = / dx sin ip (x) x 



E 



7T 

l£ 



rp2 ■ 2 <P 0) 

1„ sm — - — 



(31) 



which should be substituted into Eq. <|25|l instead of 5 SC . 

Now we turn to the case of weakly reflecting barriers, in 
which case the effective action is defined by Eqs. (|22I23|1 . 




FIG. 1: Third cumulant of the current operator S3 at zero 
temperature as a function of frequencies lu\ and u>2- Note 
that £3(0)1, U2) changes its sign with increasing frequencies 
and flattens off at |cji|, |oJ2 1 > eV. 



Combining this action with Eq. (|25() for the odd (to 
21 + 1 > 3) current cumulants we obtain 



e m TZ 

S m = — 6{ui + ...+ L0 m ) 



2tt 



2eV+ £ E 

Ml,2=±l Vj=±l 

(-I)-*- ( e y + £ ^2 coth — + Ei# 

3 = 1 7 V 3 = 1 



n 

3=1 



M1+M2 M1~M2_ 

4 H 4 2T 



(32) 



Analogously for the even (to = 21 > 4) we find 



e m TZ 

S m = 7, S(U)\ + . . . + W m ) 

Z7T 



E EM) 

Ml,2=±l Vj=±l 



2 E] w j coth 

3=1 



2T 



(33) 



x ye 

m 

*n 

3=1 



iV + > coth h > — 



3 = 1 
Hi + M2 



2T ^ 4T 

3=1 



4 ' J V 4 " 2T / 

The delta-functions <5(cji + . . . + w m ) in Eqs. 



(|29I30I32I33|1 illustrate the cumulant invariance with re- 
spect to the time shifts tj — ► tj + At. One can also 
consider "on-shell" cumulants S rn symbolically defined 



as 



S m {u>x, . . . , tjj m _i) — 

_ S m (uJi + . . . + LO m _i, — U>i, —LO2, ■ 

8{wi + . . . + u> m ) 



-(•34) 



This definition implies that one should first remove the 
S— function from Eqs. (|29I3UI32I33|) and then replace the 
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FIG. 2: Fifth current cumulant <5>s at T = 0. Two frequencies 
are fixed, Wi = u>2 = —5eV, while two others, u-j, and u>i are 
varied. 



frequencies as follows u)% — > X^JLi w 2 

The cumulants 5m, denned in 



-ux, oj 3 

—UJ 2 , ... U) m —> —U) m -X' 

this way, depend on m — 1 frequencies and in the limit 
u>i — > reduce to the standard FCS expressions for the 
current cumulants. 

The complete analytical expression for the third cumu- 
lant Ss(uji,u! 2 ) at arbitrary reflection probabilities has 
been previously derived in Ref. 8 by means of a some- 
what different approach. In the limit of small 1Z the 
result 8 reduces to that derived here. 

The third current cumulant defined by Eqs. I|32l34|) 
with m = 3 is plotted in Fig. 1 in the limit T = 0. We 
note that the overall shape of the function S^(uji, lu 2 ) re- 
mains the same at any transmission values T n (cf., e.g., 
Fig. 2 of Ref. l2l|) except for the tunnel junction limit 
(all T n <C 1) in which case no frequency dispersion of S is 
observed. In a general case the cumulant £3 becomes dis- 
persionless only at sufficiently high frequencies (or at low 
voltages) reaching a universal valu o 8 i 21 S3 — > f3e A gV/2ir. 

Higher cumulants depend on larger number of frequen- 
cies and their behavior turns out to be more complicated. 
As an example in Fig. 2 we plot the fifth current cumu- 
lant (defined by Eqs. I|32I34|) with m = 5) at T = as 
a function of two frequencies keeping two others fixed. 
Unlike the third cumulant, 1S5 does not tend to any uni- 
versal frequency-independent value in the high frequency 
limit. 



V. LOGARITHMIC CORRECTIONS AND RG 

Let us turn on electron-electron interactions. In this 
section we assume the external circuit to be purely 
Ohmic, Zs{uj) — Rs- Further we adopt the standard 
set of approximations and assume that either the bar- 
rier dimensionless conductance g — 2j^ T„, or that of 
the external circuit gs = 2n/e 2 Rs is large, j» 1, or 
gs S> 1. In this case in a wide range of parameters the 
effect of electron-electron interactions is to produce a neg- 



ative correction to the system conductance. The latter 
interaction correction is proportional to l/(g + gs) and 
depends logarithmically 1,22 on voltage and temperature, 
thus becoming large at sufficiently small eV and T. This 
divergence demonstrates insufficiency of the first order 
perturbation theory in l/(<? + gs) at low energies and 
makes it necessary to evaluate the higher order terms. 

For the sake of definiteness below in this section we 
will assume T <C eV and define 



L = In 



1 



eVR C 



where R = RR S /(R + Rs) and R = n/e 2 T n is 
the Landauer resistance of the barrier. The effective 
expansion parameter of the perturbation theory is then 
L/(g + gs)- 

Let us first evaluate the second order contribution to 
the interaction correction oc L 2 /{g + gs) 2 ■ For this pur- 
pose it is sufficient to expand the effective action to the 
third order in tp~ . One finds 2 



iS = iSs 



dt' cp~ 



1 

R 



2e 2 R 



dtidt 2 a(ti - t 2 )tp {h)<p (t 2 ) 



x {1-/3 + £cos |> + (ti)-¥> + (ia)]} 
" ft dr( V -(r))^+(r) 



6e 2 R 

27U7 
_ 3e 2 i? Jo 
x/(t2,ti)/(t3,t 2 )/(t 1 ,t 3 ), 



dTxdT 2 dT Z Lp (Tx)<p {r 2 )tp (r 3 ) 



where 



and 



12n Tn 



En 



f(r 2 ,rx) = 



Tsin [{ V +{t 2 )- V +{t 1 ))/2\ 
sinh [itT(t 2 — tx)] 

The current is expressed via the path integral 



I(t) 



JV<p 



± e iS[<p±] 



(35) 



(36) 



(37) 



(38) 



This formula can be derived, e.g., from Eq. (j3J) making 
use of the definition for the current operator 



I(t) 



e d(N L (t) - N R (t)) 
2 dt 



(39) 



where N L>R (t) = £ CT J L R , d 3 r #t ( t> r)* ff (i, r) is the to- 
tal number of electrons on the left (right) side of the 
barrier. 
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Evaluating this integral with the approximate action 
and assuming eVR C <lwe get 



I = — 



1 



O 



2I3L 

a + as 

L 3 

(a + as) 3 



2L 2 



[9 + as) 



0- 



P 2 Rs 
R + Rs 



27 



(40) 



In order to find the higher order terms of the pertur- 
bation theory in L/(g + gs) one needs to retain the con- 
tributions to the effective action of order ~ (<y9~) 4 and 
higher. An alternative way is to treat the effective action 
by means of the RG approach 3 ' 4 which allows to recover 
the leading logarithmic contributions to all orders. 

Previously, the RG equations for this problem were 
formulated either in the limit 3 gs S> 1,<? or in the oppo- 
site limit- 4 gs — > and g 3> 1. In a general case the RG 
equations for the channel transmissions read: 



dT n = T n (l - T n ) 

dL E k f k + (gs/2y 



(41) 



The implicit solution for these equations can be obtained 
in the form 



f n (L) 



T n {l-z) 



1 - T n z 

L = -^ln(l-T fe z)-^ln(l-z), (42) 



where the parameter z changes from to 1. 

Eqs. (|42|l demonstrate that as the voltage decreases 
the channels are "turned off" by interactions one by one 
depending on their transmission values. Most transpar- 
ent channels R n <C 1 remain open down to lowest volt- 
ages. Resolving Eqs. <|42l) one can explicitly determine 
the renormalized (voltage dependent) transmissions T n 
and, substituting T n into the Landauer formula, derive 
the expression for the I — V curve. Unfortunately it re- 
mains unclear whether this approach is sufficient down 
to the lowest energies/voltages in which case instanton 
effect a 6 i 23 need to be taken into account. The analysis of 
this problem is beyond the scope of the present paper. 



VI. CURRENT IN THE LIMIT TZ « 1 

Let us now evaluate the interaction correction to the 
current for a somewhat different physical limit. In con- 
trast to the previous section, here we will make no as- 
sumptions about both dimensionless conductances g and 
gs, i.e. the conductances of both the barrier and the ex- 
ternal circuit can no longer be large. At the same time 
we focus our attention on almost transparent scatterers 
assuming TZ <C 1. 

We again make use of Eq. I|38[l combining it with the 
effective action (1221 . In the path integral we perform a 
shift <p + -> ip+ + eVt, where V = V x /(l + e 2 N ch R s /Tr) is 



the voltage drop at the barrier in the absence of interac- 
tions. This shift helps to eliminate the linear in V x term 
from the action but simultaneously introduces the 
voltage V in the last term of Eq. 12/il) . The variational 
derivative of the action takes the form: 



2ieR 



(N ch -K cos (p-(t)) 



dxa(t 



tp (x) (fi (t) 
x) sin — - — cos — - — 



X cos 



eV(t~x) + v+(t)-p+(y) 



dz 



Tsm\mT(t — y)ip (z) 



sinh7rT(t — z) smhirT(y — z) 



(43) 



Here . . . stands for the terms which give no contribution 
to the current. 

Evaluating the path integral (|38J) we can put TZ = 
in the action S'[<p ± ] in Eq. I|38|) an d keep the terms oc TZ 
only in 5S sc /6ip~ . It is possible because the terms coming 
from first order in TZ correction to the action vanish. The 
integral over ip + gives S— function which fixes (p~ . Since 
the latter phase turns out to be small, \p~\ < 7r, we are 
allowed to use the action in the form 121|) . The path 
integral becomes Gaussian and can be evaluated exactly. 
The result reads 



I=-{N ch 

TT 

2eTZ 



TV 



-TZ)V 

dt a{t) e F ?{t) sin 



K(t) 



Here we have defined the functions 



K(t) = 
F(t) = 



duj 



ie 2 e~ lut 



2tt 



4loC - 



e 2 N ch 



+ 



djj 



2T 



2?r -iuC + 



sm[eVt}. (44) 

(45) 
(46) 



s^)) 



Zsi 
1 — COS LOt 



Let us point out that the form of Eq. (|44|) resembles to 
a certain extent that of the result for the I — V curve 
derived perturbatively in T n <ti 1 for externally shunted 
tunnel barriers within the so-called P(£?)-theoryi£>. Here, 
in contrast, the interaction term in Eq. l|4^|l was derived 
perturbatively in TZ — J^ n (l — T n ). Another important 
feature of our result is the presence of the term e +F ^ un- 
der the time integral in Eq. (|44p. This exponent becomes 
large in the long time limit and should be contrasted 
with the decaying exponent e~ F ^ in the corresponding 
expression^ derived in the limit T n -C 1 . 

In the limit V — > oo the I — V dependence (|44|l tends 
to the following simple form 



- iV ch 

TT V 



s e e 



(47) 



i.e. the I — V curve has the offset. This result formally 
holds for any Zs(uj) and at any temperature. In practice, 
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the offset might be difficult to observe at high conduc- 
tances and temperature. 

Below we will consider an important limit of purely 
Ohmic external environment Zs(u>) — Rs- If both tem- 
perature and voltage are sufficiently small, T, eV <C 
(R + Rs)/RRsC, the integral in Eq. (jSJl can be evalu- 
ated analytically. We obtain 



At the same time the above results are non-perturbative 
in both l/g and 1/gs and, hence, provide complementary 
information to the results obtained from the RG analysis 
discussed in the previous section. 



VII. CURRENT NOISE 



e 2 N ch 



V 



,0o 



{R + R s )\ s+g s 



2irTRR s C 
AireTIZ 



\ a+gs ) \g+gs 2nT 



sinh 



2T 



cosh ^ — cos — — 
* g+gs 



(48) 



where T(x) is the gamma-function and 70 — 0.577 is the 
Euler constant. At low temperatures, T <C eV, Eq. I|48|) 
yields the differential conductance 



d£ = e 2 N ch 

dV 7T 



Trr 1 



e 2 K f ef°(R+R s ) 
e\V\RR s C 



g+g* 



(49) 

while in the zero bias limit eV <C T one recovers the 
linear conductance 



e 2 N ch e 2 K r i 1 ~ jtgZ ) { e*> (R + R s )\^~° 



G 



2^- 



\ 2 g+g, J 



irTRRsC 



(50) 

We note that Eqs. Q44l50fl . being perturbative in TZ 1 
become inapplicable at voltages/temperatures below the 
energy scale E* , which can be estimated as 



E 



, R + Rs ( n 



RR S C \N ch 



(51) 



We now turn to the current noise and again consider 
the system with almost perfectly transmitting channels 
assuming The current noise spectrum is defined 

in a standard manner: 



5/(w) = / dtS(t)e 



(52) 



where 



S(t) = (I(t)/(0) + /(O)J(t)) - 2<I(0)) 2 . (53) 
In the absence of interactions and in the limit 1Z <C 1 
the low frequency noise spectrum scales &ai& 5/(0) oc 1Z. 
Below we will demonstrate that electron-electron inter- 
actions yield an additional contribution to 5/(0) which 
is also proportional to 1Z but can be much larger than 
the non-interacting result. In other words, a dramatic 
increase of the current shot noise by electron-electron in- 
teractions is expected provided both conductances g and 
gs are not very large. 

The noise spectrum will be evaluated with the aid of 
the path integral 



-2ie 2 / dte 



6<p- (t)6tp- (0) c 



j Vlp ± e iS[ v i 



.(54) 



This expression follows directly from Eqs. (|5I39|1 . 

The variational derivative of the effective action S sc 
can be evaluated in a straightforward manner. We obtain 



-2ie' 



6 2 S SC & ± ] 



2ie 2 K 



6tp-(t)6ip-(0) 



<p+Qt) Oa[ip-(p)]S(t) 



2e !^J± a (i) - ^Jl a (t)(coB[<p-(t)] + cos[^-(0)]) 



2e 2 TZ 

7T 



2e 2 TZ 



a(t) cos 



<p + (t) - v? + (0) - i / dz(Tcoth[wT(t - z)]+Tcoth[nTz])Lp-(z) 



Aie 2 1Z 



dx a(t — x) sin 



¥>"(*) y~(Q) 

cos — - — cos — - — 

2 2 

ip + (t) -ip + (x)-i I dz(Tcoth[nT(t - z)] + Tcoth[7rT(z - x)])<p~(z) 



: (Tcoth[7rTt] - Tcoth[7rTa;]) cos sin f-J^L 



4ie 2 K f 

/ dx a(—x) sin 

7f J 



<p + (x) - ip + (0) - i / dz (T coth{nT (x - z)}+Tcoth{nTz])Lp-{z) 



x (T coth[7rT(x — t)] — T coth[7rTt]) cos sin 



(55) 



Here we again omit terms which contribution to the path integral vanishes. 
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Since we are going to evaluate linear in TZ contributions 
to the noise spectrum, it suffices to set TZ = in the 
expression for S in the exponent of Eq. I|54|l. As in 
the previous section integrating first over tp + one can 
verify that only the values \ip~ \ < tt give a non- vanishing 
contribution. Under this approximation the path integral 
becomes Gaussian and we find 

s nt) = + ~ aS ^> 



2e 2 TZ 

IT 

Ae 2 TZ 



a(t) e 



F(t) 



COS 



K(t) 



cos[eVt] 
K{x) 



cos[eVx] 



dxa(x) e F ^ sin 
x (T coth[7rT(x — t)]+T coth[7rT(x + 1)]) . (56) 
Let us define the function 



m = F(t) + -K(\t\)sigrd 

2 I did (1 — cos coth i sin cot 
e / — : - — , (57) 



2tt 



and express the noise power spectrum as follows: 

,4<d 0„2 



Si (a) = 

2e 2 TZ 
+ 



e^TZ 2e*,„ , to 

— + (iV ch - 2TZ)u coth — 

7TO 7T Zl 



dx a(x) e*^- 1 cos[eVx] cosux 

dx a{x) e®( x > cosfeVa:] sin uix. (58) 



Aie 2 TZ coth ^ 



Making use of the property $>(t — i/T) = $(— t) one can 
prove the following identity: 



dxa(x) e*^ coswx = — — + wcoth 
20 



+ i coth — 
2T 



2T 

dx a(x) e*^ sin ws. 



(59) 



Then, combining Eqs. (JSSJl and we obtain 

2e 2 (7V ch -2^) w 
= ejeoth— + 

7T 2i 

H 2^(w±eU) coth- 



± 



2 51 



coth 



2T 

CO 

2T 



SI 



where 81 (V) = I(V) 



coth 



LU±eV 



2T 



(60) 



\N ch -K) 



V is the interaction cor- 



rection to the current, I(V) is defined in Eq. (|44|l . The 
frequency dependence of the noise spectrum is illustrated 
in Fig. 3. One can see that the noise is enhanced due to 
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FIG. 3: Frequency dependence of the noise power spectrum 
(solid line) at T = and in the presence of interactions. For 
comparison the noise spectrum in the absence of interactions 
is shown by the dashed line. Both curves are plotted for 
iVch = 1, K = 0.1, g s = 2 and V = e/2C. 



interaction at low frequencies |o>| < eV, and reduced at 
M > eV. 

It is easy to see that Eq. (|60|l satisfies the fluctuation- 
dissipation theorem (FDT). Indeed, applying an ac bias, 
V(f) — Vq coswi, and repeating the analysis of the previ- 
ous section one can verify that the real part of the zero- 
bias conductance is related to the current in the following 
way: 



Re Glui) = -I 



(61) 



Making use of this identity and setting V = in Eq. 
(|60|l . one finds 



Si(w, V = 0) = 2 [Re G(lo)] cu coth — 



(62) 



in agreement with FDT. 

In the non-interacting limit Eq. (|60|l reduces to the 
Khlus formulaic expanded in TZ to the first order. We 
also note that the result for the noise spectrum for highly 
transmitting barriers (|60|1 is to some extent similar to 
that for the case of tunnel junctions. The latter can 
be expressed in terms of the I — V dependence in the 
following way: 



SrW = eTl(-±eV\ coth""" 1 



± 



2T 



(63) 



This is a general result which holds in the non-interacting 
limits as well as in the presence of an arbitrary external 
impedance 2 ^. We observe that in both Eqs. (|60|l and 
(|63|l the effect of electron-electron interactions is fully 
described by the interaction term contained in the I — V 
curve. 
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FIG. 4: The zero frequency shot noise power spectrum 5/(0) 
as a function of voltage V at T = for N c h = 1, R — %/e 2 
and 7£ = 0.1. The conductance gs effectively controls the 
interaction strength, i.e. at small gs the interaction is strong, 
while it tends to zero in the limit gs — » oo. The interaction- 
induced excess noise is clearly observed even at rather large 
values of gs- 



Consider now the low frequency limit to — ► 0. From the 
above results we obtain 

5/(0) = (^{N ch -2Tl)+2^8I{V) S jAT 

+ 2ei^—V~SI(V)jcoth—. (64) 
In the shot noise regime T <C eV this formula reduces to 



5/(0) = 2e 



e 2 U 



V - SI(V) 



(65) 



Since the interaction correction is negative, SI(V) < 0, 
we arrive at the conclusion (cf. Ref. ,2) that for 
highly transmitting barriers the shot noise is enhanced 
by electron-electron interactions. At high voltages the 
noise spectrum takes the form 



5/(0) 



2e ir R 



V 



2c) ' 



(66) 



i.e. in this limit interactions induce voltage-independent 
excess noise which magnitude scales linearly with the 
Coulomb gap e/2C. 

In the limit eVRR s C/(R + R s ) < 1 from Eqs. flEJ 
and H48|l one obtains the result 



5/(0) 



2e 2 n(R + R s )e^ 
- r ( 2 - Fk) RR sC 



e\V\RR s C \ 1 
R + R s 



(.67) 



This expression is valid provided e\V\ > E* , where E* is 
defined in Eq. (J5IJ. In the limit g + g s — > oo Eq. ljf77|) 
reduces to the non-interacting result. For g + gs S> 1 the 



interaction effects remain weak and Eq. I|67(l agrees with 
the perturbative results 2 . In the limit g + gs 2> 1 the 
same expressions can also be reproduced with the aid 
of the RG analysis I|41I42|) combined with the formula 

Si = ^-eV J2 n Tn(l — T n ). For smaller conductance 
values the excess noise becomes larger and for g + gs of 
order one the shot noise is strongly enhanced by electron- 
electron interactions. The effect of interactions on the 
shot noise spectrum 5/(0) of highly transmitting coherent 
scatterers is clearly observed in Fig. 4. 



VIII. QUANTUM DOTS 

So far we have considered a single scatterer embedded 
in the electromagnetic environment with the impedance 
Zs(lo). Another important physical situation is that of 
so-called quantum dots. A quantum dot can be modeled 
by a system of two scatterers connected in series via a 
sufficiently small island. In this section we will demon- 
strate that some of the results derived above for a single 
scatterer can be directly generalized to the case of quan- 
tum dots. 

Let us denote the number of channels in the left and 
right barriers Ni and N r . As before, we will also assume 
that all these channels are almost open and Hi,H r <C 1. 
Throughout our analysis we will stick to a simplified de- 
scription which amounts to treating the effect of the left 
barrier on the right one (and vice versa) as that of an 
Ohmic resistor. This approximation is justified if the 
electron distribution function inside the dot remains in 
equilibrium, i.e. equal to the Fermi function with the 
temperature of the leads. This is the case either in the 
linear in voltage regime or, else, provided inelastic re- 
laxation of electrons inside the dot is sufficiently strong. 
The latter condition may apply for large quantum dots 
and/or at high enough temperatures. In addition, below 
we will ignore the gate modulation of the current flowing 
across the quantum dot. 

Let us denote the voltage drop across the left/right 
barrier as V/, r . Obviously, Vi + V r = V, where V is the 
total bias voltage. From the current conservation condi- 
tions at both barriers one finds 



/ = 



_ e 2 (N l -TZ l ) 
e 2 (N r - Tl r ) 



Vi 



v; 



-sii(Vi), 

SI r (V r ). 



(68) 



Here 5Ii is the interaction correction to the current in 
the left barrier given by the last term in Eq. 144|) and 
Eqs. H45I46|I with the following replacements: N c h — » JVj, 
l/Z s (w) -> {£)N T , C C s , where C E is the total ca- 
pacitance of the quantum dot. The interaction correction 
SI r is defined analogously. 

Resolving Eqs. (|S5|> and keeping the first order in TZi >r , 
one finds 

e 2 Nt N r T , e 2 N 2 Kl + N?1Z r , 



N 



-V 



N 2 



-V 



12 



N *XT { Nr v\-L Nl XT 



Ni_ 

N 



-V 



(69) 



where N = Ni + N r . Assuming now T -C eV and 
eVRiR r Cx/(Ri + R r ) < 1 we obtain 



it N 



2e^/ N eE c (N^ 1/N TZ l + N^ 1/N K r ) 
7r 2 iV 2-2/w r ( 2 - ^) 



7reU 
2~E^ 



1-1/jV 

(70) 



and 



dl_ 
~d~V 

e 



e 2 N t N r 



N 



2 ivr i/jv % 



N?- 1,N n r 



N 2-2/N T (2 _ 



7T 



2e'°E c 
ireV 



1/N 



,(71) 



where Ec — e 2 /2Cs is the charging energy. In the op- 
posite limit eV «T« (Ri + i? r )/i?;i? r Cs we arrive at 
the expression for the linear conductance 



G = 



e 2 NiN r 
~V N 



2 2 T(l-^)(N 2 TZ l+ N 2 n r ) (2e'*E c y / " 



r(§ 



J_) M2-1/N 
2N ) JV 



7T 2 T 



In the limit Ni = N r = 1 this result reduces to that de- 
rived in Ref. [l7| by means of the bosonization technique. 
Eq. H72|) is also consistent with the corresponding result 
of Ref. Il8l in the case of arbitrary number of channels 
Ni, r . 

Under the assumption of strong inelastic relaxation 
inside the dot current fluctuations across two barriers 
can be considered uncorrelated and one can also evalu- 
ate the noise spectrum of the system. Combining Eqs. 
SI = SVt/Rt + b, 51 = dV r /Rr +£ r and 5V t +SV r = we 
arrive at the relation between the current fluctuations in 
the system SI and fluctuations across each of the barriers 



SI 



Ni+N r ' 



(73) 



Accordingly, the noise power spectrum Si oc (SI 2 ) takes 
the form 



N 2 N 2 



(74) 



where the noise spectra Si. r arc defined by Eq. I|67|) with 
the corresponding parameters. Here we consider only the 
shot noise regime T <C eV and assume eV <C NEc, in 
which case one finds 



Si 



Ae^ (N 2 ^ + N 2 U r ) 2 (ireVy- 1/N 
tt 2 N 2 -V n T(2- jj) 6 C \2E C ) 



(75) 



As before, we observe strong enhancement of the shot 
noise by electron-electron interactions provided the num- 
ber of channels in the dot N is not large. 



Further work is needed in order to relax the assumption 
about strong inelastic relaxation inside the dot. In the 
metallic limit g ^> 1 the corresponding analysis has been 
developed in Ref. 
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APPENDIX A: TRANSFORMATION OF THE 
ACTION 

Let us demonstrate the equivalence of the actions Q 
and Q. First we evaluate the Green function Gy which 
satisfies the equation 

Gy\t 1 ,r 1 )G v (t 1 ,t2;r u r 2 ) = Sih^Sin^), (Al) 

where G v l {t,r) is defined in Eq. The general solu- 
tion of Eq. IjAlfl can be written in the form 



G v = 



-i6(ti-t 2 )Ui{t x ,t 2 ) 




+ i 





*0(t2-tl)W 2 (tl,i2) 

pv{t) -pv(t) 



Wi(ti,*) 

U 2 (h,t) J V Pv(t)-1 -Pv{t) 

Ux{t,t 2 ) 

o u 2 {t,t 2 ) 



(A2) 



Here pv(t) is an arbitrary operator. Below we will 
use the convention according to which the product of 
operators involves only the coordinate integration, i.e.: 
U 2 {t u t 2 ,ri,r 2 ) = JdVW(ti,f2,ri,r / )W(ti ) t2 ) T- , ,r 2 ). 
It is also important to keep in mind that < t\, t 2 < t, 
and the Green function Gy implicitly depends on the 
final time t. 

The task at hand is to evaluate the operator pv(t). For 
this purpose let us make use of the Dyson equation 

Gv(ti,t 2 ) — G (ti,t 2 ) 



dt'G (ti,t') eV(t') G v (t', t 2 ), (A3) 



where 



V(t) 



Vx{t) 
V 2 (t) 



(A4) 



and Go is the Keldysh-Green function of noninteracting 
electrons: 

r _ ( -iO(ti-t 2 )U (t u t 2 ) \ 

' 1 i0(ia-ti)Z4>(*i,*2) J 
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W o (*i,0) 
Wo(*i,0) 

«o(0,t 2 ) 
U (0,t 2 ) 



Po -po 
po — 1 -po 



(A5) 



Here Uq is the evolution operator for non-interacting elec- 
trons defined by eq. JSJ with Vj. j2 = and po is the initial 
reduced single-electron density matrix operator defined 

as (r\po\r') = tr ^^(r')4f^(r)poj , where po is the ini- 
tial many-particle density matrix. Substituting the gen- 
eral solution (|A2I) into Eq. I|A3|I and making use of the 
Dyson equation for the evolution operators 

Wi, 2 (M') = U (t,t')+iJ dsU (t,s)eVi, 2 (s)Ui, 2 (s,t'), 

we find2&: 

pv(t)=Ui{t,0)po 

x [l + (W 2 (0,t)Wi(t,0)-l)A>r 1 W2(0,t). (A6) 

As a next step, let us fulfill the following replacements 
eVi -> eF+ + AV-/2, eV 2 -> eV+ - AF"/2 and evaluate 
the derivative idS e i/dX, where iS e i = 2TrlnGy 1 . With 
the aid of Eq. (@J) we obtain 



J ds J d 3 r(G^ A ( S , S ,r,r) 
-G 2 v 2 x (s,s,r,r))eV~( S ,r). (A7) 



,dS el 



Employing Eqs. i|A2[l. (|A6I) and using the properties of 
the trace of the product of operators, we find 

dS i 



+ U 2 ,x(s,t)p v , x (t)U 2 , x (t, S ))eV-(s 



i I ds 

x tr 



(Z/ 2i a(0, t)U hX (t, s)eV-(s)U 1 ,x(s, 0) 

+ W 2 , A (0, 8)eV-(a)U2,\(«> *)Wi,a(«, 0)) 
x p [l + (W 2 , A (0,i)^ A (i,0) - l)p ] _1 



(A8) 



What remains is to integrate Eq. (| A8|) over A from to 
1. This task is accomplished with the aid of the identity 

^■(U 2 ,x(0,t)U 1 ,x(t,0)) = 



ds(U 2 ,x(0, t)U hX (t, s)eV-{s)U h x(s, 0) 

+ W 2 ,a(0, s)ey- (s)U 2 ,x(s, t)U!,x(t, 0)) . (A9) 

Since the action S equals to zero at A = 0, we arrive at 
the final result 

2TrlnG^ 1 = 2trln[l + (JA 2 (0, t)U x (*, 0) - l)p ] ,( A1 0) 

which proves the equivalence of Eqs. iJSJ and Q). 



APPENDIX B: FULLY TRANSMITTING 
BARRIERS 

Let us demonstrate that the action So (|18|1 can be ex- 
actly transformed to the form l|19fl . For this purpose we 
first expand the action i|18f) in the following series 



iS [a, (f + ] = 2Tr< a(x)p (y-x)e 



(-I)V 



7V= 2 



iV 



(Bl) 



where 



Jn = j dxi . . .dxN po(xi - x 2 )a(x 2 )po(x 2 - x 3 ) . . . 

x p (x N -i - x N )a(x N )p (x N - xi)a(x N ), (B2) 



and rewrite Jn in the form 



Jn = I dx I — — 



dE 
2^ 



du>\ . . . duj n 



\N 



x n(E + u>{)n(E + uj± + w 2 ) . . . n(E + k>i 
x a(w!)a(cj 2 ) . . . a{u N ) e i (<"i+-+ w *)* 



n{E) 



+ Wjv_i) 

(B3) 



where a(u/) = J dxe iux a(x) is the Fourier transform of 
a(x), and we have also introduced the cutoff at a large 
negative energy -E min . 

We will use the following property of the Fermi func- 
tion 



n(E)n(E + oj) = - (l 



coth^) n(E) 



<-uth—)n(E + u). (B4) 



This equation is analogous to the Ward identity for 
the electron Matsubara Green function: Q(uj + f2,p + 
q)g(u,p)=g(n,q)[g(u,p)-g(u + Q,p+q)]. The latter 
identity ensures that all higher order symmetrized po- 
larization bubbles for the one-dimensional electron gas 
vanish thus turning the RPA approximation into an ex- 
act procedure^. Analogously, the identity 1B4I) helps to 
simplify the integrals (|B3ll . 

In what follows we will employ the procedure similar to 
that applied within the imaginary time formalism^!. Let 
us interchange {u>i, ■ ■ ■ ,u>n} N — 1 times in the following 
way. We first put lon in front of lo\, then we place it 
in-between lji and u> 2 , and so on. This set of changes is 
summarized in the table: 



{uji UJ 2 W 3 W4 ... LJ N } 

I 

(1) {U N Ul UJ 2 UJ 3 ... LO N -l} 

(2) {UJI WJV W2 w 3 ... UJ N -l} 

(3) {uji U> 2 UJ N ^3 • ■ ■ U N -x} 



(B5) 



(N-l) {ui uj 2 



lun- 2 lon Wff-i}. 
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Having made these changes in the integral (|B3|) we ex- 
press Jn as a sum of N — 1 corresponding integrals di- 
vided by TV — 1. Afterwards we apply the identity l)B4jl in 
each of these integrals excluding lon from the arguments 
of the Fermi functions. Specifically, in the term corre- 
sponding to the sequence (1) in the table l|B5|) we split the 
product of the first two Fermi functions n(E)n(E + Wjv), 
in the term (2) we split the product of the second and 
the third Fermi functions, and so on. Then we get 



Jn 



1 



N - 1 



, dE f du>i . . . dwjv 



-Emin 



2tt 



(27T) 



TV 



x [f(-u N )n(E) + f(w N )n(E + cj n )} 

x n(E + uj n + u)\) . . . n(E + u N + Wi + . . . + Wjv-2) 

x 5(wi)5(«a) . . . d(uj N ) e i (^+-+»^ x 

dE f du>% . . . dux 

- > i i/.r i 

N - 



i y. 1 i r°° dE r 



(27T) 



N 



n{E)n(E + w\) . . . x [f{-u N )n{E + + ... + cj fc _i) 

f(uj N )n(E + uji + ... + Uk-i + u! N )} 

x n(E + u)i + ... + uJk-i + Wjv + Wfc) x . . . 

X 7l(E + CJi + . . . + UJk-l + Wjv + Wfe + ■ • • + Wjv-2) 

x 5(wi)a(w2) . . . 5(wjv) e i( - u ' 1+ - +u '^ x , (B6) 

where we have defined f{uj) — \ (l + coth 

Now let us perform the frequency shifts u>j — ► + w^v 
in all the terms as well as the energy shift 75 — > E + lun in 
the second term of the first integral in order to eliminate 
u>n from the arguments of the Fermi functions. As a 
result we arrive at the following expression 



J 



1 



N 



N — 1 



dx 



dE 



divi 



, dux 



x [/(-wjv)a(wi - wjv) + /> J v)&(^i)e MNX ] 
x n(E)n(E + uji) . . . n(E + uj x + . . . + uj N - 2 ) 

xd(uj 2 )---&(uJN)e liu ' 1+ - +u ' N - l)x 

1 *z} 



N 



k=2 



dE 



(27T) 



x [f(-WN)a(wi,-i)a(wk -&n) 
+ f(uj N )d(u) k -i - w N )a(uj k )] 
n(E)n{E + wx)..- n(E + u x + 



+ 0JN-2) 



x a(wi)a(w 2 ) . . . a(w fc _ 2 )a(^fc+i) • ■ ■ a(wjv) 



a i(aii + ...+w N _i)x 



1 



TV - 1 



dx 



--En 



2^ 



\N 



x n(E)n(E + bJi) . . . n(E + uj x + . . . + uj N - 2 ) 
x f(u N )a(vi)afa) . . . a(u N ) e **+~+>>«)*. (B7) 

We will assume that |wi|, . . . , \u>n\ "C E m [ n which implies 
that n(E)n(E + u>x) . . . n(E + u>i + . . . + LUN-2) = 1 in 
the interval —E m i n < E < —E m i n + lun- This fact allows 



us to evaluate the integral over E in the last term. In 
addition, one observes that /(wjy) + I{—^n) = 1 and, 
hence, the function /„ should drop out from all the other 
terms. We obtain 

I f , f°° dE f du>\ . . . cLun 



Jn 



dx 



N 



x d{u>\) . . . cl{ujm) e* 



(27T) 
N-l , 



3=1 



x n(E)n(E + u) X ) ■ ■ ■ n(E + w\ + . . . + ui N ^ 2 ) 



1 



dx 



du)\ . . . duj 



N - 1 

x a{u)x)a{oj 2 ) ■ ■ ■ a(uj N ) e i ( w i+-+ w ")* (B8) 

Returning to the coordinate representation and assum- 
ing a(x) — * for |.t| — * 00, we find 

1 



Jn 



dx\ . . . dxN-i po(xi 



x 2 



N-l 

x Po{x 2 - x 3 ) . . . p {x N ^ 2 - x N ^ 1 )po(x N ^ 1 - xx) 

N-l / N-l 



j'=i 



JV-l, 



x)a(y) + a(x)a N 1 (y) 
4(7V- 1) 



dxdy . .a 



To summarize, making use of Eq. (|B4|) we have re- 
duced the number of the functions pa under the inte- 
gral by one. One can employ this procedure further and, 
again applying Eq. i|B4|) . finally reduce the number of 
the functions po down to two. As a result, Jn takes the 
form 

. N-l 
J N = dxdyp (x - y)p (y -x)^] ~j kN a k {x)a N ~ k {y) 
J k=l 

+ f ^a(x-y)f^(3 kN a k (x)a N - k (y).(B10) 

J I 1 



Now let us apply the identity 

Po(x - y)po(y -x) = p o (0)S(x -y)- v — y> , (Bll) 

which follows directly from Eq. l|B4p . and regroup the 
terms in Eq. (|B10|) . Then we obtain 

Jjv = J dxp o (0)u N a N (x) 

+ ( ^ a{x - y) Y^ akNa k {x)a N- k{y) _ (B12) 

J 71 k=l 

In order to find the coefficients un and a k N let us 
rewrite Eq. (|B9|) in the form 

d JN-^aixy^^} 



Jn[ci(x)] = -i 
dxdy 



d\ N-l 

a N - 1 {x)a(y)+a(x)a N - 1 {y) 



2n 



a(x - y) 



4(JV- 1) 
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and derive the following equations: 
N- 2 1 



OCN-1.N 



N - 1 
N -2 



OL\ jv-1 



4(JV- 1) : 



fe-i 



N-l-k 



= Wiv-i (Bl 

together with the boundary conditions iti = 1 and ai^ 
— 1/2. These ennations ea,n be solved with the result 



'Mi — 1) OtkN 



N 



~ Ak(N - k)' 



(B14) 



Let us now recall that the function po (y — x) in Eq. 1|18|) 

is multiplied by e l 2 . Hence, in Eq. (|B12|> po(0) 
has to be replaced by 

lim/5 (y-x)e 4 ^ =po (o) + _. 

I-^S 47T 

This observation brings us to the following expression for 
Jn- 

Jn= J dx L(0) + 

Substituting J/v into Eq. (|B1J) . we finally arrive at Eq. 

APPENDIX C: WEAK REFLECTION LIMIT 

Let us derive the effective action of a coherent scat- 
terer in the weak reflection limit making use of the per- 
turbation theory in R n . In order to proceed we define the 
operator 



D a (x,y) = (x\ [1 + p a] po\y), 



(CI) 



where (x\p \y) = p (x-y) and (x\pa\y) = p (x-y)a(y). 
Expanding Eq. i|13|) in f''f' one finds 

iS { J ] = -2K f dx 



lim I D ai (x,y)e l * M ) ai(x) 

x^y 



+ lim D a2 (x, y)e 



8ft J dxdy D a2 {x,y)D ai (y 7 x) e 



a 2 (x) 

i[ip+(x)-ip+(y)] 



. <p (x) . U3 (y) 
x sin Y „ v 1 sin Y } yj . 



(C2) 



Here we have defined oi )2 (i) = 8(t-x)8{x)(e ±lv {x) -l), 
and D a (x,y) = (x\D a \y). 

The function D a (x,y) can be found in the following 
way. According to Eq. (|18[) one has 



iSo[a + da, 0] = iSo[a, 0] + 2 j dx D(x,x)8a(x) 

dxdy D(x,y)5a(y)D(y,x)5a{x) + 0{Sa 3 (xj). (C3) 

On the other hand, from Eq. I|19fl we find 

iSo[a + 8a, 0] = iSo[a, 0] — f dx Po ^\ Sa 2 (x) 

J (l + a{x)y 

+ I dx 2p ^ + ^^^ X - y)Hl + am Sa(x) 



dxdy 
An ' 



',(x - y) 



1 + a(x) 

Sa(x)Sa(y) 
(1 + a{x))(l + a(y)) 



Sa 2 (x) ln[l + a(y)} Sa 2 {y) ln[l + a(x)} 



2{l + a{x)f 2(1 + a(y)) 2 

+ 0(5a 3 (x)). (C4) 

Comparing Eqs. (|C3|I and l|C4(l we obtain 



D a (x,x) 



po(0) + ±f dya{x - y) ln[l + a{y)] 



1 + a(x) 



(C5) 



D a (x,y)D a (y,x) 



Po(0) 



dy 

—a(x-y) ln[l + a(y)] 

4-7T 



a(a; - y) 



47r(l + a(a!))(l + o(i/)) 



(l + a(x)) 2 
(C6) 



These relations are consistent with the following form of 
D a (x,y): 



D a {x 7 y) 



■ J dz(coth7rT(:r-z)-cotli7rT(t/-z)) ln[l+a(z)] 



y/(l + a(x))(l+a(y)) 
x p (x-y) . 



(C7) 



Since Eqs. (|C5IC6(1 do not uniquely determine the func- 
tion D a (x,y), Eq. I|C7|I has to be additionally checked. 
To this end we make a shift a{x) — > a(x) + Sa(x) and 
find the linear in 5a(x) correction to D a (x,y). From Eq. 
(|C1|) we obtain 



5D a (x,y) 



dzD a {x,z)8a{z)D a {z,y). (C8) 



Defining the function f(x) = \5{x) — |Tcoth7rTx, and 
making use of the identity 

Po(x-y)po(y-z) = [f(x-y) + f(y-z)]p (x-z), (C9) 



which is a direct consequence of Eq. |B4J| , one can verify 
that the property (|C8|) also holds for the function l|C7(l . 
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We then conclude that the two functions IjClfl and (|C7|) 
may differ only by a shift of the argument a(x). The 
latter shift is zero since in both cases at a(x) = one 
gets D a= o(x, y) — pa(x — y). Thus, we conclude that Eq. 
(IC7|) is indeed correct. 

What remains is to substitute this result into Eq. (|C2|) 
and get 



+ 81Z dxdyp (x - y)p (y - x ) e ^ + ^)-v + (y)] 



(CIO) 



o 

x e J* Q * dz[T cothnT (x-z)-T cothnT (y ~ z)]<p~ (z) 



iSg* = -411 / dxp (0)[l-ca&<p-(x)] 
Jo 

■A / ax(p (X) sin cp (X) 

* Jo 

HP* 

H / dxdya(x — y)ip (x)sunp (y) 



In the last term of this equation we again apply the iden- 
tity IjBlljl and arrive at the final expression for the effec- 
tive action (J2HJ. 
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